Interplay between the Fulde-Ferrell phase and Larkin-Ovchinnikov phase in the 
superconducting ring pierced by an Aharonov-Bohm flux 
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We study the phase diagram of a superconducting ring threaded by an Aharonov-Bohm flux 
and an in-plane magnetic Zeeman field. The simultaneous presence of both the external flux and 
the in-plane magnetic field leads to the competition between the Fulde-Ferrell (FF) phase and 
the Larkin-Ovchinnikov (LO) phase. Using the Bogoliubov-de Gennes equation, we investigate the 
spacial profile of the order parameter. Both the FF phase and the LO phase are found to exist stably 
in this system. The phase boundary is determined by comparing the free energy. The distortion of 
the phase diagrams due to the mesoscopic effect is also studied. 
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I. INTRODUCTION 



In recent years, one of the inhomogeneous su- 
perconducting states, known as Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state, has received a lot of inter- 
est. This superconducting state with periodical spacial 
variation of order parameter (OP) was first proposed in- 
dependently by Fulde and Ferrell^ and by Larkin and 
Ovchinniko\s2, in 1960s. The possible evidence of its ex- 
istence has been reported in certain unconventional su- 
perconductors^ and the possibility of its realization in 
trapped cold atomsj^iii^'^ In literature, the state is col- 
lectively known as the FFLO state.— Actually, they are 
two kinds of states with slight difference: the order pa- 
rameter of the LO state is real and spatially inhomoge- 
neous, which breaks the translational symmetry, while 
the order parameter of the FF state has a uniform mag- 
nitude, but an inhomogeneous phase similar to that of a 
plane wave, breaks the time-reversal symmetry. Accord- 
ing to previous studies, the FF state is usually unstable, 
and unfavorable in comparison with the LO state. Al- 
though the LO to FF phase transition was predicted in 
Ref. \m a more recent studyi^ shows that there is no 
stable FF phase in such a system and there is no LO to 
FF phase transition cither. The authors in Ref. [l^, men- 
tion a possible FF state in a momentum space study, but 
as to the best of our knowledge, a realization of stable 
FF state in the presence of a Zeeman field has not been 
reported yet in a real space calculation. 

As is well known when a Zeeman field is added to a 
superconductor, the LO state becomes favorable in com- 
parison with the BCS state, irrespective of the geometry 
of the superconductor. Meanwhile, we notice that in a 
superconducting ring, which is threaded by a magnetic 
flux, the Aharonov-Bohm (AB) flux breaks the time re- 
versal symmetry in much the same spirit as that in the FF 
phase As a result the FF state comes out. An interest- 
ing question is then if we add both the magnetic flux and 
an in-plane magnetic field, how will the two phases com- 
pete with each other? Motivated by this observation, we 
study in this paper the interplay between this AB flux- 



driven FF phase and the Zeeman field-induced LO phase. 
It is of great interest to study the phase transitions and 
phase diagram in such a system. The investigation is car- 
ried out in a tight-binding model for a superconducting 
ring pierced by an AB magnetic flux, and in the presence 
of a Zeeman magnetic field. We solve self-consistently 
the Bogliubov de Gennes equation for the superconduct- 
ing order parameter and determine the phase diagram by 
comparing the total energy. We find that for this system, 
there are four different phases when we vary the two pa- 
rameters, magnetic flux <I> and the Zeeman field h. More 
interestingly, we also study the mescscopic effect. 

The paper is organized as follows: in Sec. Ull we intro- 
duce the tight-binding model and present the mean-field 
treatment. In Sec. IIII[ we numerically carry out the cal- 
culation of superconducting order parameter as a func- 
tion of the magnetic flux and Zeeman field, and deter- 
mine the phase diagram by comparing the free energies. 
Section HVl is the discussion and conclusion. 



II. MODEL AND MEAN-FIELD TREATMENT 

We consider a one-dimensional superconducting ring 
threaded by an external magnetic flux $ (see Fig. [1]). 
Meanwhile, there is an in-plane magnetic field B, which 
generates the Zeeman spliting and gives rise to the in- 
homogeneous pairing. The system is described by the 
following Hamiltonian 
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Here Uj — iy c'^'^*/^*", where is the bare hopping 
coefficient and $o = hc/e is the normal-state flux quan- 
tum, and N is number of lattice sites for the ring; c]^ 
(cicr) is the creation (annihilation) operator on the i-th 
lattice site with spin cr = ±1 for spin up and down elec- 
trons, arising from the interaction between the magnetic 
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is solved by iteration. Here T is the temperature the 
Boltzmann constant feg = 1 has been taken.) . Notice 
that the quasiparticle energy is measured with respect to 
the chemical potential. 




FIG. 1: (Color online) Schematic illustration of the setup. 
A superconducting ring is threaded by an external magnetic 
flux, denoted by $. A magnetic field B is applied in the plane 
of the ring. The ring is connected to the ground to ensure 
that the chemical potential is fix, but the electron number 
may fluctuate 



field and the spin of the electrons; rii.o- 



is the 



particle number on the i-th site with spin cr; g is equal 
to 2; is the Bohr magneton and B is the strength 
of the in-plane magnetic field; V is the strength of the 
on-site pairing interaction; /i is the chemical potential. 
For simplicity, we define h = g^sB as the strength of 
the Zeeman field. In the present work, we take tij to be 
t between nearest neighboring sites and zero otherwise. 
Within the mean-field approximation, the Hamiltonian 
([T|) is reduced to 



vY^iA4^ci + h.c.)+Y: 
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III. NUMERICAL RESULTS 

In our numerical calculation, we take the energy unit 
t = 1, and the chemical potential fi = —0.5, the interac- 
tion strength V = 2, and the ring size = 50. Though 
the system size is far from the thermodynamic limit, it 
already gives the same phase boundary as that of infinite 
A^. The order parameter structure depends not only on 
the Zeeman field h, but also the magnetic flux We 
notei^ii^ that all physical quantities have already been 
a function of 4> with a period of $o even in the normal 
state. Therefore, it is sufficient for us to consider the 
magnetic fiux in the range € [0, <f>o]. In the absence of 
the magnetic flux $ = 0, the Bardcen-Cooper-Schrieffer 
(BCS) order parameter A = 0.351 for h = 0, and the 
LO state is stable for hd < h < hc2 with hci = 0.23 
and hc2 = 1.56. The system becomes normal (A = 0) 
for h > hc2- In the presence of the magnetic flux, the 
magnetic flux can induce a change in the structure of 
the BCS state in an s-wave superconductor, namely a 
crossover from the BCS state in the absence of a mag- 
netic flux to a FF state with a magnetic flux when the 
Zeeman field is low h < hd- When the Zeeman field 
increases, the LO becomes favorable and both BCS and 
FF states give in. If we continue to increase the Zeeman 
field, the amplitude of the pairing potential of the LO 
phase will be suppressed by the Zeeman field until it dis- 
appears finally, and the system enters the normal state. 
In the following, we will numerically construct the phase 
diagrams. 



where Ai = V {ci^Cii) is the pair potential. To diagonal- 
ize this Hamiltonian, we employ the following Bogoliubov 
transformation 
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corresponding to the eigenvalues where 7j/ and 7^ are 
the quasi-particle operators. The coefficients (uf^,, wf^.) 
satisfy the Bogoliubov-de Gennes (BdG) equationii^ 
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where Hi 



-tij — fi6ij + crhSij, and Hi 



potential 



- ah*Sij. The self consistent equation of the pair 
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A. Phase boundary in plane 

We first focus on the low temperature case (3 = 1/T = 
200 (corresponding to T = 0.005). In the absence of the 
magnetic flux, there are three different phases: BCS, LO, 
and normal. In the presence of the magnetic flux, there 
are also three phases, FF, LO, and normal state. In the 
following, we study the phase transitions and the phase 
boundaries for fixed temperature when varying h and 

In order to check if the FF state becomes the ground 
state, we assign a periodic phase to the order param- 
eter at each site as an initial condition. Similarly, we 
assign a constant phase to see if BCS state becomes the 
groimd state. For a set of fixed parameters (/i, TV, T), 
different stable solutions (with different order parameter 
textures) could be obtained from different initial config- 
urations. For example, one may find both stable LO- 
typc OP and FF-typc OP for the same set of parameters 
(A^ = 50, /3 = 200, $ = 0.25$o, h ^ 0.25). Even there 
are more than one stable LO type solutions for the same 
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set of parameters, which means different net momentum 
of the Cooper pair. To distinguish one state from other 
competing states (including BCS state and FF state), we 
choose the energetically most favored one by comparing 
their free energies. For the model H]), the free energy is 
given by 



F 
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h) 
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Here, we just compare the summation of the first two 
terms, because the third term is a constant for all solu- 
tions of different phases. In the following we determine 
the phase boundary between FF state, LO state, and 
normal state. 



1. First order transition between the FF and LO phases 

When determining the pair potential self-consistently 
by iteration, we find that in certain range of the strength 
of the in-plane Zeeman field h, different initial configu- 
rations of the pair potential lead to different stable so- 
lutions. In another word, there are more than one sta- 
ble solutions through iteration. For example, when we 
fix (f> = $o/4, and vary the magnetic field in the range 
0.08 < h < 0.29, stable solutions of both the FF type and 
the LO type can be arrived at through iteration. The free 
energies of these two types of stable solutions are listed 
in Table I. It can be seen that the LO state becomes en- 
ergetically favorable when the magnetic field is equal to 
or greater than hci = 0.21. In addition, the free energy 
at hci is continuous, but its first order derivative is not 
continuous. Hence, we conclude that for a fixed magnetic 
flux 4> = $o/4, there is a first-order phase transition be- 
tween the FF and LO states at hd. Similarly, we fix 
magnetic flux $ at different values and we can find the 
threshold value of h at which the system changes from 
the FF state to the LO state or vise versa. Thus for a 
fixed temperature f3 = 200 and fixed system size = 50, 
the phase transition line between FF and LO state is de- 
termined by comparing the free energy of the FF phase 
and the LO phase, and we plot it in Fig. [21 To ensure 
that the phase boundary given by = 50 is close to that 
of the thermodynamic limit, we change the system size 
to = 200, and we find the phase boundary does not 
change. For iV = 50 and N = 200, the magnitude of the 
OP Ai in BCS phase is the same. Hence the result based 
on = 50 can be regarded as in thermodynamic limit. 
It can be seen that the first-order transition line is not 
parallel to the $ axis, so we can turn the flux to make 
the system change from the LO phase to the FF phase 
or vise versa. We call this phase transition AB effect in- 
duced phase transition. We can also see that the phase 
boundary between the LO and FF states is symmetric 
around <!> = $0/4, and the period of FF phase is $0/2. 



TABLE I: Free energies (up to a constant — X^^i(/i-f^)) for 
stable solutions of FF state and LO state. Here the ring size 
is TV = 50, the magnetic flux $ = $o/4, and the temperature 
(3 — 200. It can be seen that there is a first-order phase 
transition from the FF state to LO state when the in-plane 
magnetic field is tune across h = 0.21 





h = 0.20 


h = 0.21 


h = 0.22 


h = 0.23 


h = 0.24 


h = 0.25 


FF 


-56.1487 


-55.6487 


-55.1487 


-54.6487 


-54.1487 


-53.6487 


LO 


-56.1019 


-55.6319 


-55.1619 


-54.6919 


-54.2219 


-53.7519 
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FIG. 2: (Color online) Phase diagram of the superconduct- 
ing ring in the plane. Here the ring size is A'^ = 50, and 
the temperature is T = 0.005 (/3 = 200). Notice that the 
boundary line between the FF and LO phases has the peri- 
odicity in "I> with a period of $0/2 while that between the 
LO and normal state phases has the periodicity in $ with a 
period of $o- Specifically, the LO state exists in the range 
[0.23, 1.56], [0.21, 1.52], and [0.23, 1.51] for $ = 0, $o/4, and 
$0/2, respectively. 



2. Second-order transition between the LO and normal 
state phases 

If we continue to increase the in-plane magnetic field 
above the value h = 0.29 for $ fixed at $o/4, all initial 
configurations of the pair potential will lead to the LO 
state, or only the LO state becomes stable. Meanwhile 
the amplitude of the pair potential decreases and the pe- 
riod of the modulation of the pair potential is shortened 
continuously. Further increase of the Zeeman field leads 
to the reduction of the pair potential until it vanishes 
gradually. When the magnetic field reaches hc2 — 1.52, 
the amplitude of the pairing potential vanishes, or the 
LO state is completely depressed by the in-plane mag- 
netic field, and the system changes from the LO state to 
the normal state. If we do the iteration from zero pair 
potential, we will find that when h > hc2 the stable solu- 
tion for the pair potential is zero (normal state). When 
h ^ hc2, the stable solution is an LO state. There is no 
coexistence area of the LO and the normal states in the 
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h axis. Hence we conclude that the phase transition at 
hc2 for a fixed 4> is a second-order phase transitoin. Our 
result is consistent with previous studies.— i^i^^ 



B. Phase boundary in h^T plane 
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FIG. 3: (Color online) Phase diagram of the superconducting 
ring in the h-T plane. Here the ring size is A'^ = 50, and the 
magnetic flux are $ = (a), and $ = $o/4 (b), respectively. 
For $ = 0, the LO state emerges at low temperature T < 
0.12 and relatively high magnetic field 0.23 < h < f.56. For 
$ = '&o/4, the LO state emerges at low temperature T < 0.11 
and relatively high magnetic field 0.2f < h < 1.52. The LO 
state to normal state transition (Black with open circles) is 
of second order, while the BCS (or FF) to the normal state 
transition (Red with open squares) is of first order. The zero- 
field transition temperature is around Tc = 0.21 and 0.20 for 
$ = and $0 /4. 

In the preceding subsection, we study the phase tran- 
sitions when we vary the magnetic flux $ or the in-plane 
magnetic field h. The temperature is fixed at a very low 
value. Hence these phase transitions can be regarded as 
quantum phase transitions. In this subsection, we will 
study the phase transitions induced by thermal fluctu- 
ations, and determine their phase boundaries. We will 



fix the magnetic flux $ and vary the temperature /? or 
the in-plane magnetic field h. First we consider the 
case in the absence of the magnetic flux $ = 0. We 
fix the temperature at /3=20 (T=0.05), /3=10 (T=0.10), 
/3=6.67 (T=0.15), [3=5 (r=0.20), and /3=4 (T=0.25) re- 
spectively, and do the iteration separately. The phase 
boundary between the BCS and LO states is determined 
in a similar way to that in Sec. IIIII A. It can be seen that 
when we fix the magnetic flux to be zero, and tune the in- 
plane magnetic field or the temperature, the system will 
change between the BCS, LO, and normal states. As can 
be seen from Fig.[3I^a). the LO phase emerges below the 
critical temperature T « 0.12. We note that the BCS to 
the LO state is first order, and the LO to normal is sec- 
ond order. When the magnetic flux is nonzero, the BCS 
state will be replaced by FF state with the phase dia- 
gram, as shown in Fig.[3jb), very similar to the zero-flux 
case. 



C. Mesoscopic effect 

Another interesting question is the mesoscopic effect. 
In this subsection we will study the mesoscopic effect 
by fixing the temperature and decreasing the ring size. 
As mentioned in the above discussion, the ring size 

= 50 already gives the same phase boundary as that of 
A'^ — > oo. A simple check is that when we increase the ring 
size to A^ = 100, and A^ = 200, we find that the phase 
boundaries do not change in comparison with that for 
A = 50. This means that for the current model, A^ = 50 
can be treated as in the thermodynamic limit. However, 
if we decrease the ring size, for example to A = 20, the 
mesoscopic effect will occur. First, in the plane, the 
LO phase will shrink dramatically and the FF phase will 
expand (see Fig.[4lja)). This is because (1) the influence 
of the magnetic flux on the system will increase and the 
influence of the in-plane Zeeman field will decrease rela- 
tively, and (2) having a finite size restricts the periodicity 
of the LO order parameter. At a given Zeeman field, if 
the period is not commensurate with the corresponding 
ring, solutions of the LO state will have to be modified to 
be commensurate with system size, which results in some 
energy cost. Therefore, the LO state will shrink. Second, 
the periodicity of the magnetic flux changes from <E>o/2 
to $0- This is because the system size is so small that 
the Cooper pair can no longer be treated as a whole, and 
can only be treated as two separate electrons. We show 
in Fig. m^b) the phase boundary for A^ = 10 and the 
re-entrant behavior of various phases can be seen in the 
phase diagram. 



IV. DISCUSSION AND CONCLUSION 

Based upon a tight-binding model, we study a one- 
dimensional s-wave superconducting ring subject to an 
in-plane Zeeman field and a magnetic flux by solving the 
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FIG. 4: Phase diagram in the plane. AU the parameters 
are the same as that in Fig. [5] except that the ring size is 
= 20 (a) and N = 10 (b). The empty area represents the 
normal state. The gray area represents the FF state, and the 
area covered by the thin black lines represents the LO state. 
The thick black lines represent the BCS state. It can be seen 
that with the decrease of the ring size, the FF phase expand 
a lot and the LO phase shrink dramatically. The period of 
FF phase also changes from $o/2 to $o 



In the absence of the magnetic flux, there is no FF phase, 
and the Zeeman field induces the transitions between the 
BCS, LO, and normal states, which has been studied ex- 
tensively,ilii^ii2iSai2ii2^i^i^i2^i^ Our results agree well 
with the previous studies in a two-dimensional system 
that the energetically favorable state for s-wave super- 
conductor is a one-dimensional stripe-like LO state. This 
suggests the first-order transition between the BCS and 
LO states while a second-order transition between the 
LO to normal states. Our study goes beyond that and 
indicates a stable FF state due to the magnetic flux. The 
mesoscopic effects are also studied. When the system size 
decreases, two mesoscopic effects arise: (1) the LO phase 
in the h-^ plane shrinks, and the FF state expands due 
to the enhancement of the Aharanov-Bohm effect; (2) 
the periodicity of the external magnetic flux will change 
from (f>o/2 to $o. 

The following remarks are in order: (1) Though 
we study a one-dimensional model, the system should 
not be regarded as a mathematically one-dimensional. 
The current study can be easily extended to the two- 
dimensional and other geometry, such as a torus config- 
uration threaded by a magnetic flux. It can be expected 
that a similar phase transition between LO state and FF 
state will occur. (2) For the one-dimensional case, the LO 
state exists in a broader range of parameters space [h — T 
space, see Fig. [5]) than that of two-dimensional and three- 
dimensional cases fiSiSli^ which makes it easier to access 
experimentally. (3) In Ref . [2^. it is reported that an trap 
potential with arbitrary configuration can be achieved. 
Hence, we expect that the result presented in this pa- 
per should be able to be observed experimentally in cold 
Fermions under current experiment technique. (4) In the 
thermodynamic limit, N —>■ oo, the FF state reproduces 
the BCS state, because the phase gradient of the order 
parameter is vanishingly small. This result agrees with 
our intuition that when the ring size becomes infinity, 
the influence of the magnetic flux can be neglected. (5) 
The effect of the impurity is not included in the current 
study, and will be given in our future studies. 



BdG equation in real space. In the presence of a mag- 
netic flux, a crossover from the BCS state to the FF state 
is obtained when the in-plane magnetic field is not very 
strong. If we increase the strength of the in-plane mag- 
netic field, the LO state becomes favorable, and a FF to 
LO phase transition occurs. With the further increase 
of the in-plane magnetic field strength, the magnitude of 
the pair potential of the LO state is suppressed, and dis- 
appears finally with the system entering the normal state. 
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